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ABSTRACT 

We present spectral analysis from XMM-Newton observations of the Vela pulsar. We analyzed thermal emission from 
the pulsar dominating below ~ 1 keV since extracted spectra are heavily contaminated by nebular emission at higher 
energy. Featureless high-resolution spectra of the Reflection Grating Spectrometer aboard XMM-Newton suggest the 
presence of a hydrogen atmosphere, as previously indicated by Chandra results. Both the temperature and radius 
are consistent with those values deduced from Chandra. The derived Chandra and XMM-Newton temperature of 
T°° ~ (6.4-7.1) x 10 5 K at its age of ~ 10 4 years is below the standard cooling curve. 



INTRODUCTION 



The Vela pulsar is one of the best-studied pulsars in the radio band. It has an 89 millisecond period and the spin- 
down parameters indicate that the Vela pulsar is about 10 4 years old with magnetic field 3 x 10 12 G. Nevertheless, 
multi-wavelength observations in the past two decades showed the Vela pulsar is a strong source at all wavelengths 



from the radio to the gamma-ray band. In the X-ray band, thermal emission was discovered by 7?0&4r(Qgelman et al. 



1993). ROSAT observation also revealed a peculiar pulse profile with three peaks in the soft X-ray band. Non-thermal 



emission from its surrounding nebula is dominant at higher energies. 

Recent Chandra observations revealed even more complex properties of the Vela pulsar an d its nebula. With its 



high angular resolution, a Chandra image showed detailed structure in the surrounding nebula ( Helfa nd et all 2001 



|Pavlovetal.ll2001bh . No spectral features were found in the high energy resolution Chandra LETGS spectra of the 
resolved pu lsar emission . No d etection of spectral features excludes the presence of heavy element atmospheres on 
the surface ( Pavlov et al. , 2001 ah . 

The Vela pulsar is also a strong source of non-thermal emission in the optical, X-r ay and gamma-ray ba nds. The 
non-thermal component from the pulsar was resolved for the first time by Chandra dPavlov et al. , 2001ah . Recent 
RXTE observations revealed complex timing properties asso ciated with non-therma l emission from the p ulsar, show- 
ing energy-dependent multiple peaks in folded lightcurves dStrickman et al. , 19991: Harding et al. , 2002 ). However, 
the origin of the peaks in the soft X-ray band, whether thermal or non-thermal, has yet to be determined. 

The XMM-Newton observation of the Vela pulsar was performed with a 90 ksec exposure time on 1 Dec 2000. 
XMM-Newton is sensitive to faint sources s uch as isolated neutro n stars due to its large effective area. Since the 
nebula turned out to be compact ~ 10-20" (Helfan d et all 12001). the imaging capability of XMM-Newton cannot 
resolve the pulsar from the nebula. Therefore, the thermal component dominates at E < 1 keV and the nebula emission 
is dominant at E > 1 keV in XMM-Newton spectra. We could not spatially resolve the non-thermal emission from 
the pulsar or study different components of the nebula. Instead, using the XMM -Newton data, we did the following: 
(1) a spectral fit to high sensitivity Europian Photon Imaging Camera (EPIC) (Striid eret all 1200 H data was used 
to determine thermal properties with high accuracy, ( 2) a search for spectral features in high resolution Reflection 
Grating Spectrometer (RGS) Jden Herder et all 1200 lh data was perfor med. All the data were processed by version 
5.3.3 of the Scientific Analysis System (SAS) pipeline dLoiseau 1 12003). 
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Fig. 1. XMM-Newton/EP\C-PU spectra fitted by Blackbody+Power-law (BB+PL) (top) and Hydrogen 
atmosphere+Power-law model (H+PL) (bottom). Residuals at ~ 0.5 keV are due to instrumental Oxygen 
edge. 



EPIC PHASE-AVERAGED SPECTRAL ANALYSIS 

We analyzed EPIC-PN data to determine the global features of the X-ray spectrum. Source data were extracted from 
a circular region with radius r = 1' and events with CCD pattern were selected. Such a source area includes more 
than 90% of the source counts at energies up to ~ 5 keV. However, the extracted spectra are heavily contaminated by 
nebular e mission at E 1 keV . Unlike the recent Chandra results which resolved the pulsar emission from its nebula 
emission ( Pavlov et al. . 200 lab , our analysis focuses on thermal emission from the pulsar at E < 1 keV. Systematic 



errors from nebular contamination are discussed later in this section. Background spectra were reduced from an 
annulus region between r = 1' and r = 2'. Total EPIC-PN count rate from the r < 1' region after subtracting 
background was 27.0 cts/sec. 

We rebinned the data so that each spectral bin contained more than 25 counts. For the following spectral analysis, 
we used ready-made Response Matrix Files (RMFs) and Auxiliary Response Files (ARFs) generated by the SAS tool 
"arfgen". For spectral fitting, we adopted a two-component model comprised of a thermal and non-thermal (power- 
law) spectrum. We applied two different thermal models to t he data: a bla ckbody model (BB+PL model) and a 
magnetized hydrogen atmosphere model provided by V. Zavlin (Za vlinLl2003h at B = 10 12 G (H+PL model). When 
we fitted a hydrogen atmosphere model, we fixed the neutron star mass and radius to 1.4M Q and 10 km. Both thermal 
models fit the data well, yielding \ 2 — LI (Figure 1 and Table 1). Figure 2 shows an Nh-T°° contour plot with fixed 
power-law index and radius. Hydrogen atmospheres are in general harder than a blackbody since free-free absorption 
dominates in the X-ray band. Therefore, the temperature fitted by a magnetized hydrogen atmosphere model is lower 
than a blackbody fit (hence a larger radius for the hydrogen atmosphere fit). 

The derived distance (d = 256 pc) is consistent with the distance from the parallax measurement (d = 294^g[j 
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Fig. 2. Contour plot in N H -T°° plane. Blackbody (top) and hydrogen atmosphere at B = 10 12 G (bottom). 



pc) ( Caraveo et all | 200lh. For hydrog en atmosphere models, T°° and R°°/d depend on other input parameters (e.g. 
M and R) weakly (Za vlin et al.Lll998l) . The apparent radius corresponding to the measured distance is inferred to be 
-R°°/^294 = 15.0^5^ km where £^94 is the distance in units of 294 pc. 

We note that non-thermal spectra from the pulsar (obtained by Chandra observation) are about two orders of 
magnitude weaker than the nebular spectra, so the power-law component observed by XMM-Newton is dominated by 
the nebular contribution. The fitted power-law index (7 ~ 1.6) probably represents a superposition of different nebular 
components, which were resolved by the Chandra observation ( Pavlov et al. . 2001bh . The estimated luminosity of the 
power-law component in t he EPIC -PN data is consistent with the Chandra result for the nebula (L ne {, = 8.3 x 10 32 



ergs s _1 for d = 294 pc) dPavlov et alll2001bh . 



The derived parameters from both the blackbody and magnetized hydrogen atmosphere fit are in good agreement 
with the Chandra results. However, since the angular resolution of XMM-Newton is poorer than that of Chandra, 
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Table 1. Fitted parameters to XMM-Newton EPIC-PN 
spectra (pulsar plus nebula) in comparison with Chan- 
dra data (pulsar only). 

Blackbody + Power-law 
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N H [10 20 cm" 2 ; 
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Magnetized H atmosphere c + Power-law 
EPIC-PN Chandra 
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Errors for EPIC-PN data are 1 sigma level including both 
statistical and systematic uncertainty. For luminosity of non- 
thermal component in the range of 2-10 keV (L p i), the error 
refers only to the statistical uncertainty since the smaller ex- 
tracted region picks up less nebular emission. 
a We assumed the distance is 294 p c from the latest parallax 
measurement dCaraveo et all 1200 lb. We resc aled the Chan- 
dra results in Table 1 in Pa vlov et all (12001 ah to d = 294 pc. 
b Note that these parameters are different from the Chandra 
results in comparison. However, the fitted parameters to the 
XMM -Newton data is domi nated by the nebular component, 
while Pavlov et al. extracted spectra from the pul- 



sar. 

c We fixed the neutron star mass to 1.4M Q and the radius to 
10 km for spectral fitting. Here B = 10 12 G. 



our fitting parameters can be affected by systematic errors due the the nebular contamination. These errors are, in 
general, larger than the statistical ones. In order to estimate such errors, we examined the variations in the hydrogen 
atmosphere and blackbody parameters by using as extraction regions annuli of various radii; the errors reported in 
Table 1 reflect both statistical and systematic uncertainties. We can see that, when considering the effect on the nebular 
contamination, Chandra and XMM-Newton temperatures are in agreement at 1 -sigma level. 

SEARCH FOR SPECTRAL FEATURES IN RGS DATA 

Reflection Grating Spectrometer data were taken in "spectroscopy" mode in order to take full advantage of the high 
spectral resolution. First, we extracted the 95% Point Spread Function (PSF) of the source region for both RGS 1 
and 2 data. The background spectrum was reduced from outside the 98% PSF region. We generated responses by the 
SAS command "rgsrmfgen". Total count rate was 1.0 cts/sec for both RGS 1 and 2 data. We extracted smaller regions 
of the source (50% and 70% PSF region) to maximize the number of thermal photons compared to nebular emission 
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Fig. 3. XMM-Newton/RGS spectra fit by Hydrog en + Power-law model. The spectra were extracted from the 
70% Point Spread Function (PSF) region. 

photons. Figure 3 shows RGS spectra extracted from the 70% PSF region and fitted by the H+PL model. The data are 
plotted in the RGS bandpass (A = 6-38 A) after rebinning wavelength grids to A ~ 0.1 A. 

The spectral fit by both the BB+PL and the H+PL model yield xl - 1-2 for both RGS 1 and 2 data. The fitted 
parameters were similar to the EPIC-PN results in Table 1. Several deviations from the continuum fit with less than 2 
sigma significance are seen, but none of these wiggles were found in both RGS 1 and 2 data simultaneously. Therefore, 
we conclude that the RGS spectra are featureless. 

SUMMARY 

The results obtain ed with XMM-Newton for the Vela pulsar are in agreement with previous findings by Chandra 
dPavlov et all 1200 lallbl) . The RGS data indicate that the thermal spectrum is featureless, and the EPIC-PN spectrum 
is well fitted by a magnetized hydrogen atmosphere model with B = 10 12 G. As it has been found with Chandra, 
with respect to a fit with the simple blackbody, the atmospheric model has the advantage that it gives a radius for 
the emitting region compatible with neutron star equations of state. Also, when Chandra data are fitted with an 
atmospheric model, the index of the power-law component above ~ 1 keV is remarkably consistent with that of 
the PL observed at optical and hard X-ray wavelengths. The power-law component observed with XMM-Newton is 
dominated by the nebular emission and the value of L p i ~ 8 x 10 32 erg s _1 is again consistent with that observed by 
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Chandra JPavlov et allbOQlbh . 



The presence of hydrogen is not sur prising since a tiny amount of hydrogen constitutes an optically-thick layer in 
the atmosphere dZavlin & Pavlovu2002h. Hydrogen could be accumulated on the surface by accretion or generated by 
spallation of fallback material (IjBildsten et al. , 1992h . The derived large radius from the H+PL model implies that the 



thermal emission originates from the whole surface. 

There is an important implication for the neutron star physics which is consistent with the Chandra and XMM- 
Newton results. Surface temperature for a given age reflects cooling processes via neutrino emission in the core. 
Recent studies show that neutron star cooling curves are dependent on various factors such as proton-to-neutron 
ratio, neutron star mass, magnetic field strength, neutron superfluidity and presence of exotic matter. Among various 
proposed models, a model which assumes normal composition (n, p and e) and the indirect URCA process was often 
quoted as the standard cooling model. The standard cooling model fits most cooling neutron stars, althou gh recent X 



ray ob servations found several neutron stars have surface temperatures below the standard cooling curve (Slane et al 



2002h . For the Vela pulsar, the derived Chandra a nd XMM-Newton te mperature of T°° ~ (6.4-7.1) x 10 5 K at its age 



10 4 years is below the standard cooling curve ( Tsuruta et al. . 2002h . 
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